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Abstract 

We present thermodynamics of the quantum su{l, 1) Landau-Lifshitz model, foUowing our earher 
exposition [J. Math. Phys. 50, 103518 (2009)] of the quantum integrability of the theory, which is based 
on construction of self-adjoint extensions, leading to a regularized quantum Ifamiltonian for an arbitrary 
n-particle sector. Starting from general discontinuity properties of the functions used to construct the 
self-adjoint extensions, we derive the thermodynamic Bethe Ansatz equations. We show that due to non- 
symmetric and singular kernel, the self-consistency implies that only negative chemical potential values 
are allowed, which leads to the conclusion that, unlike its su{2) counterpart, the su{l, 1) LL theory at 
T = has no instabilities. 
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1 Introduction 



There has been a renewed interest in the Landau-Lifshitz (LL) model in recent years, mainly due to its 
intcgrability, and, thus, its role as a toy-model to check the gauge/string duality as both having underlying 
integrable structures (see for recent reviews Oil])- A large amount of literature has been devoted to this 
subject, mostly concerning its classical integrability aspects. At this point, there is relatively complete 
understanding of the classical integrable properties of the LL model [3]-[5]. The quantum integrability, on 
the other hand, has been shown to be surprisingly complex and non-trivialQ Firstly, as demonstrated by 
Sklyanin |7], one has to clearly distinguish the su(l,l) and the su{2) models, as they have completely 
different physical properties. This is due to construction of the positive-definite metric in the corresponding 
Hilbert spaces. So far, only the sit(l, 1) case has been addressed in literature from this point of view, as in 
the 5^(2) case the construction of the positive-definite metric seems to be a quite complicated mathematical 
task [8] . The next interesting issue arises in the anisotropic LL model, where the algebra of observables has 
to be modified for consistency to be a quadratic one. While natural from the lattice regularization point of 
view [QllTO] , it is not clear how to derive it directly within the framework of the continuous integrable model. 
And finally, the LL model is distinctly different from other known integrable models due to its highly singular 
nature in the quantum mechanical description Indeed, it gives rise to S"{x) types of interactions which 
are not mathematically easy to deal with |T2]. As a result, the quantum Hamiltonian was not possible to 
obtain until recently. 

In our previous paper [B], we have solved one of the outlined above problems, and by constructing the 
self-adjoint extensions with the help of the unusual vector-like states in Hilbert space, specified by the non- 
trivial scalar product, we had explicitly presented the quantum Hamiltonian, and derived the spectrum of 
the su(l, 1) isotropic model, in complete agreement with all known results. In the current paper we work 
out the thermodynamics of the model, relying on our previous construction. We re-derive the Bethe Ansatz 
equations, find the integral equation describing the excitation and finally construct the thermodynamics, 
following the standard methods. We show that unlike other models, in this case one has to start from a 
finite temperature case, and only then consider a careful T — > limit to avoid contradictory results. One 
interesting consequence of the thermodynamic Bethe Ansatz equations for the sm(1, 1) LL model is that 
there is no Fermi point, which means that only the negative values for chemical potential are allowed. We 
use this feature to argue that the 5-^(1, 1) field theory at T = has no instabilities, unlike the perturbative 
results [13] for the su{2) model. 

We emphasize that we use LL model as a testing ground for such singular potentials, as the above 
mentioned aspects and difficulties of the LL model can be equally transferred to other interesting models, 
e.g. the fermionic AAF model [14j . quantum mechanical description of which has the same highly singular 
S"{x) behavior. Our goal here is to learn to deal with quantum integrability directly in the continuous 
case, without resorting to the lattice regularization schemes, which can be a daunting task by itself. We also 
consider only the sw(l, 1) isotropic LL model, leaving the essential complications of the su{2) and anisotropic 
cases to a future work. 

Our paper is organized as follows. In section 2, we give a short description of the Landau-Lifshitz model, 
and give only essential details of our method of constructing the self-adjoint extensions and give the quantum 
Hamiltonian, as well as the exact spectrum. On this basis, in section 3, we use the discontinuity properties 
of the functions used to construct the Hilbert space, re-derive in this way the Bethe equations, and following 
this initial setup work out the thermodynamics. Finally, in conclusion, we remark that the results of the 
current paper open several interesting questions, and we outline possible directions. 



^ Sec an extensive list of literature in [6] 
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2 The LL model, quantum Hamiltonian and self-adjoint exten- 
sions. 

The classical Hamiltonian of the full anisotropic Landau-Lifshitz modelH has the form: 

e 



H 
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(5,S,5,S)+472((53)'-l)] (1) 



Here the scalar product (S, S) = (S'^) -e (S"!) -e (S^) = 1 for the vector S = (S'^, 5*2, S^) depends on the 
parameter e = ±1, which gives rise to two different models corresponding to the su(l, f ) and su{2) cases. In 
what follows we will consider only the su(l, 1) (e = 1) and isotropic (7 = 0) case, as there are, as we noted 
above, principal difficulties with su(2) or anisotropic cases in the quantum theory. The classical dynamics 
is governed by the Poisson structure: 

{S^(x),S^{y)}^±tS^{x)6{x-y) 

(2) 

{S-{x),S+{y)} = 2tS'{x)S{x-y) 
which is replaced, in the quantum theory, by the commutation relations for the S'-operatorsl^ 

[S^(x),S^(y)]^±S^{x)6{x-y) 



(3) 



[S-{x),S+iy)] ^2S^{x)Six-y) 



The representations of ([3]), in the ferromagnetic vacuum S^{x)\0) = |0); S (x)|0) = 0, can be constructed 
using the vectors of the form 

\fn) = / dxi...dXnfn{xi..Xn)S~^{xi)...S^{Xn)\0) (4) 



where /n(xi...x„) are continuous and decreasing sufficiently fast functions for the integral (|4]) to be well 
defined. It can be shown that only the su{l, 1) case guarantees a positive definite scalar product (g„|/„). 

The problem one generally faces in continuous integrable models is that the quantum trace identities, 
defined via the standard procedure of first finding the transfer matrix T{u) and then taking the trace, lead, 
unlike the classical case, to polynomials which contain operator product at the same point, and therefore are 
not well-defined. This issue, however, was neglected for a long time in the majority of the classical integrable 
models, as for example the Thirring model endowed with the d{x) type of potential, by using various 
(mathematically not well defined) procedures of taking regularizations such as e{x)S{x) = 0. In the LL 
model, as well as in the fermionic AAF model [T4 , this cannot be done anymore, and formal manipulations 
with singular functions lead to meaningless expressions. As a result of these principal difficulties, the quantum 
Hamiltonian for the LL model, as well as other charges, were not possible to extract from the trace identities 
until recently. 

In our previous paper we partially solved this problem, by finding the regularized continuous quantum 
Hamiltoniar0 and giving the complete description of the resulting Hilbert space. We refer the reader to the 
publication [B] for complete details, and only present here the necessary formulas and essential results. There 
are essentially two steps involved in our construction. First, we regularize the continuous Hamiltonian by 
the split-point method so that the quantum Hamiltonian has the form: 

Hq = limiJ, (5) 



^We use the notations of [7] 

^ This is true only for 7 = 0. In general, the algebra (|3]l should be modifiod to a quadratic Sklyanin algebra [9]. 
^ We will present a method to construct the regularized charges in all orders in the upcoming paper |15| 
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where 

He = \j dudvF,{u, v) [-duS^d^S^ + duS+d^S- + dud^ {S^{u)5{u - v)) - dudj{u - v)] (6) 

Here the function F{u, v) is any sufficiently rapidly decreasing smooth function, which makes the integral 
well-defined, and the dependence on an arbitrary parameter (or set of parameters) e is not essential as long 
as we require that 

limine (w, v) = S(u — v) (7) 

This condition simply means that the quantum Hamiltonian Hq is local. It can be shown that the action 
of the Hamiltonian Hq on an arbitrary n-particle state 



|/n)= j dTx f{x) S+{x^)---S+{x^m (8) 

has the form 

HqIU) = - I d^{Af{'^))Y[S+{xm (9) 

(here the Laplacian A = J^iLi 9f) leading to the following matching conditions: 

- {d,f{x) - ddml';Zt\ = dAfi^)U=-, > > J- (10) 

This was shown to be equivalent to self-adjointness of some Hamiltonian H (for the explicit expression for 
H see [B], here we will not need it) acting on the Hilbert space generated by vectors of the type: 

/ /i(^i) \ 

f2{xi,X2) (11) 

with the scalar product, defined, e.g. for n — 2 case as 

oo 

= 1 J g*i{x)fi{x) dx + jj g;{x,y)f2{x,y) dxdy (12) 

— oo a^T^y 

which was shown to reproduce beautifully Sklyanin's earlier attempts to deal with the ill-defined object in 
the LL model. The matching conditions (jlOp were also shown to correspond to the the n-particle S -matrix 
factorization property 

Sn{Pl, ■■■,Pn) = J|S'2(k,Pj) (13) 

where the two-particle scattering S-matrix for the su(l, 1) LL model can be derived from for n — 2 and 
has the form: 

2{pi -p2)+ipiP2 

2(P1 -P2) ~ WlP2 
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Let us also note here, that in terms ol the cluster fields (xi) and ^'^(xi) (see for details [7]) the state 
(jlip for n = 2 can be written in the form 



P1P2 



(15) 



dxidx2 (c(pi,p2)e'(P^^^+P^^^^ +c(pi,p2)e^(^'^^=^+P^^^^) ^'+(xi)«'+(x2) 



2(Pl-P2)+»PlP2 



|0) 



where c{pi^p2) = ^^''\(p^^,'J)^^^ ■ Thus, the 5'-matrix has the meaning of scattering of the quanta correspond- 
ing to the ^'^(a;) cluster field. Note also, that the first term is usually absent in other models. 

To sum up, our method, being mathematically strict, gives clear systematic description and understanding 
of the properties of integrable systems endowed with singular potential. 



3 Thermodynamics 

Having obtained the spectrum and the jS-matrix, we proceed to thermodynamics of the model, following the 
standard methods (see, for example, [Tl ] ll6 [ IT7 ] ). One can easily repeat the analysis of [11], and we present 
only the final results and emphasize the differences with the non-linear Schrodinger case. Note, that while the 
analysis of the basic thermodynamic properties is quite similar to other models, there are several interesting 
features that require a careful analysis. In particular, we show below that one cannot start directly from the 
T = case, as it leads to contradictory results. Instead, due to singular behavior of the kernel, one has to 
consider the general T ^ case, and then take the T — > limit. Note also, that essential differences will 
arise when considering the correlation functions at finite temperature, as it would involve the properties of 
the correct Hilbert space ([TT|) with the scalar product p^ . 

Using the results of the previous section, we obtain the Bethe equations for the su(l, 1) LL model in the 
system of the length L : 

^^P,L ^ -rT2(K^P^2±™ (16) 

Note that the S'-matrix for the su{l, 1) that we consider here is different from the one obtained in su{2) 
by a factor "-2" in front of {p, - pj) [13 H (This is not surprising, considering these are two different 
models.) While the absolute value of this factor is not of a great importance, the sign plays a big role, 
essentially leading to different physics: it immediately follows from (fT6)) that unlike the su{2) case, in the 
su(l, 1) theory the momenta p are always real and as a consequence the theory does not have bound statesl^l 
It is convenient to introduce variables u = l/p0 This is because in the general anisotropic case, the Bethe 
equations have simpler form in terms of the u variable, which is associated with the spectral parameter in 
the inverse scattering method. Thus, (jl6l) can be also written in the equivalent form: 

L ^ 

h Y^Q{uj - Uk) = 2TTnj (17) 

Uj 



where Q{u) — zln( ^^^^^ ), and the numbers nj — n'^ + "^^r^i G J = l...A^. Since the kernel is 



K{u) EE 9„e(u) = - \ ^ (18) 
1/4 + 



^The perturbative calculations for the su{2) model lead to the S-matrix: S = FT — ^jJJPiPj. 
®It can be shown that this result is true in general for the anisotropic sn{l, 1) LL model as well. 

'^It is obvious already from that pi = is a special point in momentum space. For the moment, we will consider this 
space as R\{0}, so that the M-variables are well-defined, and extend to the full R later. 
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we see that the function Q{u) is a monotonicaUy decreasing function of u. This is contrary to the behavior 
of the closely related (in fact, equivalent on the classical level) non-linear Schrodinger model, for which the 
corresponding Q function is a monotonicaUy increasing function of momenta. Nevertheless, the solutions 
can be shown to exist and unique due to existence of the convex Yang- Yang type of action: 

Syy — ~^^^lnitj + '^'^''^^n-jUj — / 'd{v)dv (19) 

Using (|17|) and ([TBI) one can obtain the estimate 

2'n(ni — ni.) 1 1 2TT(ni — ni.) 2tt 

— — — > =Pi ~Pi> — -7^ > r 20) 

L - Uj Uk ^ ~ L{1 + D) - L{l + D) ^ ' 

This is exactly the same estimate as in the case of the NLS model, which allows us to consider the thermo- 
dynamic limit of the model. To proceed, one introduces the density of states in complete analogy as 

piPk) - 77 ^ 7 (21) 

L{Pk+i -Pk) 

and then takes the thermodynamic limit L — > cxd; N — > oo; D = N/L ~ (const). Because the zero 
temperature case is rigorously obtained as a limit of the T ^ system, we will first consider this general 
situation, i.e., the derivation of the thermodynamic Bethe Ansatz (TBA) equation. The initial steps of the 
construction are parallel to those in [11; for the NLS system, and we only mention one interesting difference, 
that in the process of the minimization of the partition function one has to be more careful, since the effective 
kernel K{pj,pk) = --(^K{pj,pk) = ^ Yj^^iillj-ilpk)^ appears in the equations in our case is not 

symmetric under the exchange of the momenta pj and p^. The resulting TBA equation has the form: 

oo 

e{p)^p^-h-^ J dfiKi^i, p) ln(l + e-^^'^)/^) (22) 

— oo 

At this point we should recall that the point p — Q has been excluded so far from our momentum space. 
Thus, we have to verify that the solution to ((22|) is continuous at this point. This will lead to a severe 
constrain on the chemical potential h. To demonstrate this, let us note that the naive conclusion that 
e{p) — > —h as p — > (which seems to follow from the fact that K{ij,,0) = 0) is not correct. In fact, we 
have the following behavior of the kernel when p — > ±0: 

\im K{f,,p)^27rS{fiT0) (23) 

Using ([23]) in ([22]), we obtain 

e(±0) = rin(e-''/^-l) (24) 



First consequence of ((24|) is that if a solution exists, then it is continuous at the point p = 0. Now wc have 
two distinct cases: 

i) h> 0. In this case a solution does not exist because (e"'*/-^ — l) is always negative or zero. 

ii) h < 0. This range of the chemical potential is allowed. So we can proceed with the same iterative 
proof as in the case of NLS to establish the existence and uniqueness of the solution of 



Here we have one of the most striking differences between LL and NLS models: while in the latter, all the 
values of h are allowed (although the thermodynamics, and especially T ^ limit, drastically depend on its 
sign), in the former, only negative values of the chemical potential are allowed. 
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Before we move to discuss the case of zero temperature, T = 0, it is instructive to compare our situation 
with NLS when parameter c (the strength of the deha-functional potential) is going to zero. This is relevant, 
because all the specifics of our model is due to effective turning- off the interaction when one of the momenta 
is approaching zero. In the case of the NLS model, the TBA equation takes exactly the same form (|22l) . 
where the kernel now is given by 

Once again, the naive setting c to zero would lead to a theory of free bosons. But it is easy to see that this 
is not the case for c — > +0, because of the property of Kc{^tP) analogous to 

lim KcU,p) = 2TrS(fi-p) (26) 

c->+0 

which allows to solve TBA exactly with the solution being 

e(p) =Tln(e(f'-'')/^-l) (27) 

(Notice the striking similarity between (p4|) and ((27|) !) The analysis of this equations also demonstrates clear 
difference of the system for /i > and h < 0. 

Now we give a brief consideration of the zero temperature limit for the case of the allowed chemical 
potential, h G (— cx),0). From (|24l) one can immediately see that in the zero temperature limit 

£(0) = -h>0 

Using the standard iterative proof of Yang and Yang (see comment after ((24| ). we easily establish that this 
is the minimum of the function s{p). As a consequence, e(p) is always positive for T — oH Then, from ()22|) 
we find that in this limit we have 

e{p) ^p^ -h>0 (28) 



In fact, Eq. (|28|) is exactly what one has in the case of NLS model for h < and T = 0, so all the results of 
the NLS model are applicable here. The main difference is that the other sector of h, namely, the h G [0, oo) 
sector is missing in the su(l, 1) LL model. 

To construct the excitations, we should again start from the general T ^ case. Using the continuous 
form of (|17|) one can write the analogue of the Lieb integral equation: 



1+ J K{p,\)pp{\)d\ = 2npt{p) (29) 

— oo 

where, as usual, Pp{p) and pt{p) are the distribution densities respectively of particles and vacancies. 

The excitations are obtained from (|17p in the standard manner, we refer the reader to the monograph 
for complete details of the similar derivation for the NLS model. For example, a one-particle excitation 
with the momentum Ap is described by the equation: 

oo 

-e(l/p-l/Ap) + 2^F(p|Ap)- j K(p,p)-lMhl}_df, = (30) 

— oo 

where the shift function F{p\Xp) = describes transition from the set of initial momenta {p} to the set 
of excited momenta {p}. Similarly, one can easily write the equations describing the hole excitation as well. 



^Let us remind, that in the su{2) case the perturbativc calculations I13| lead to instability. The main assumption there is 
the quantum integrability, which was verified in the first non-trivial order in I18| . 
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Using ([5U1) . the excitation energy takes the forn^ 



oo 



AE{Xp) = £o(Ap) - / e'M 



F(M|Ap) 



dfi 



(31) 



1 + e^t^y^ 



where the energy of the free particle eo{p) = — h,. It can be shown that AE{p) = e(p), where the latter 
function is determined from the TEA integral equation (1^^ . As usual, this demonstrates that e{p) is indeed 
the particle excitation energy. The proof consists of manipulating with TBA equation (j22l) and equation 
([50)1 . In particular, using the equation (1^51) . it is easily seen that AE{p) and e{p) have the same p^ — h 
asymptotic behavior in the T — >■ limit. 

Finally, we only mention here that the computation of the correlation functions at the finite temperature 
(rirl S'^{x)S{y) \^t) requires a separate investigation, as the unusual Hilbert space structure (fTT|) and the 
scalar product (I12|) will essentially modify the results. We leave this to a future investigation. 

4 Conclusion 

We have presented the essential steps to construct the thermodynamics for the isotropic su(l, 1) LL model, 
which is an interesting toy-model due to its singular potential structure. Wc only used our method of 
constructing the self-adjoint extensions and regularized Hamiltonian to derive the Bethe Ansatz equations, 
as well as the S'-matrix factorization property. This enables us to derive the spectrum and proceed to 
derivation of the thermodynamics. Several key features make the thermodynamics of the LL model an 
interesting system. Firstly, we showed, using non-symmetric and singular behavior of the kernel, that the 
self-consistency of the TBA equations restrict the chemical potential to be a negative function. This leads 
to a rather restrictive character of the spectrum. We also showed that the su{l, 1) theory is stable, unlike 
its su{2) counterpart. Another essential conclusion was that one has to start from the T ^ and only then 
consider a careful zero-temperature limit. 

One interesting question, which we mention just briefly, is to to compare the finite size corrections for 
T ~ for the energy with the corresponding corrections at infinite length and finite temperature. Let us 
remind, that, e.g., in the NLS model it can be shown than in the first case the finite size corrections have 
the form: 



where R = 1/T, and vp is the Fermi velocity. This is to be compared with E' = TF{T)/L , at the limit 



where in the second term Cy is the heat capacity at constant volume. The reason these two expressions do 
not coincide is, of course, the lack of Lorentz invariance in the NLS model. However, these computations can 
be easily repeated for the relativistic AAF model, and we expect these two expressions to coincide [T9_. On 
the other hand, in the case of the LL model we encounter a problem, since the Fermi point, and therefore, the 
Fermi velocity are not defined. In this case, one has to consider the finite size correction at finite temperature 
and then consider the L — > oo limit. We leave this to a future investigation. 

We again emphasize that the LL model should be considered as a toy model to understand the nuances 
of the singular potentials. It will be interesting, as a next step, to apply the results of this paper to the 
AAF model [14], which shares many general features with the LL model, but which is more interesting 
due to its relativistic invariance, as well as non-trivial bound states structure. Other interesting and related 

^It can be seen from l|30p that the function F(p\Xp) has a discontinuity at p = 0. However, it is easy to see that the integral 
I I31II is well-defined. 




(32) 





(33) 
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developments of our approach include generalization of our method to the anisotropic LL model and, what 
should be more challenging due to problems mentioned above, to the su{2) LL model. All these problems 
are currently under investigation, and we hope to report on the results in the future. 
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